Measurement strategy for spatially encoded photonic qubits 
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We propose a measurement strategy which can, probabilistically, reproduce the statistics of any 
observable for spatially encoded photonic qubits. It comprises the implementation of a two-outcome 
positive operator-valued measure followed by a detection in a fixed transverse position, making 
the displacement of the detection system unnecessary, unlike previous methods. This strategy 
generalizes a scheme recently demonstrated by one of us and co-workers restricted to measurement 
of observables with equatorial eigenvectors only. The method presented here can be implemented 
with the current technology of programmable multipixel liquid-crystal displays. In addition, it can 
be straightforwardly extended to high-dimensional qudits and may be a valuable tool in optical 
implementations of quantum information protocols with spatial qubits and qudits. 

PACS numbers: 42.50.Dv, 03.67.Hk 



I. INTRODUCTION 

Optical implementations of quantum information tech- 
nologies based on the spatial degree of freedom of sin- 
gle photons have at their disposal several approaches 
to embody a qubit or qudit structure into this contin- 
uous variable. For instance, one can encode the qubit in 
the spatial parity of the photon's one-dimensional trans- 
verse modes [1] or one can truncate (via filtering) the 
infinite-dimensional discrete Hilbert space of Laguerre- 
Gaussian [2] or Hermite-Gaussian [3] modes in which the 
photon's two-dimensional transverse distribution may be 
expanded. In simpler approaches, however, the qubits 
and qudits can be encoded either in the discrctizcd 
transverse position-momentum of the photons [4-6] or 
their angular-position-orbital-angular-momcntum vari- 
ables [7]. Here, we consider this former type of spatial 
encoding where the discretization of the one-dimensional 
transverse modes is achieved by making the photons pass 
through an aperture with D slits [e.g., see Fig. 1(a) for 
D = 2] [4, 5]. Hereafter we will refer to them as spatial 
qubits and qudits. 

Recently, considerable effort has been directed toward 
finding methods to measure spatial qubits and qudits. 
This is a fundamental issue for any possible application 
since the extraction of information from a quantum sys- 
tem requires the performance of a measurement. In this 
respect, Neves et al. [8] have shown that the spatial fil- 
tering in transverse positions at the near and far field 
of a double slit corresponds to projections in the poles 
and equator of the Bloch sphere [Fig. 1(b)], respectively. 
This enabled characterization of the entanglement [8, 9] 
and quantum tomography [10, 11] of two-photon spatial 
qubits. Later, Taguchi et al. [11] went a step further 
and showed that any vector over the Bloch sphere can be 
measured by moving a "pointlike" detector in a trans- 
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verse position between the near field and the far field of 
the double slit, both achieved with the help of a lens sys- 
tem [Fig. 1(a)]. Despite successfully accessing the whole 
Bloch sphere, this method is difficult and time consuming 
in practice since it requires moving the detection system 
for each measurement. Also, the detection efficiency de- 
creases when the detector moves away from the diffrac- 
tion envelope peak(s), which demands artificial compen- 
sation in the count rates when a pair of basis vectors 
is measured in asymmetric positions around the optical 
axis. 

A different approach for measuring only the equatorial 
states of spatial qubits and qudits has been recently de- 
veloped by one of us and co-workers [12]. There, phase 
shifts were imprinted in each slit of the array (by a spatial 
light modulator) and after this the photon was detected 
in a fixed transverse position in the focal plane of a lens, 
corresponding to an equally weighted superposition of 
the "which-slit" states without any relative phase. This 
enabled the postsclcction of noncanonical mutually un- 
biased basis (MUB) and the realization of the so-called 
MUB tomography [12, 13]. 

In this work, we generalize on this latter scheme and 
propose a measurement strategy which can, probabilis- 
tically, reproduce the statistics of any observables, not 
only those with equatorial eigenvectors. The strategy 
is as follows: First, we implement a two-outcome pos- 
itive operator- valued measure (POVM), which, with a 
given success probability, transforms the spatial qubit 
state into a convenient form. Then, the successfully 
transformed qubit is detected in a fixed transverse po- 
sition. By recording and normalizing the count rates for 
two such measurements, corresponding to two orthogonal 
projectors of [11], we can reproduce the statistics of the 
observable. Thus, the method presented here requires 
neither moving the detection system nor compensation 
in the count rates and can be implemented with the cur- 
rent technology of multipixel liquid-crystal displays. In 
addition, it can be straightforwardly extended to higher 
dimensions and may be a valuable tool in optical implc- 
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mentations of quantum information protocols with spa- 
tial qubits and qudits. 



II. SPATIAL POSTSELECTION 
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Assume that a source produces a monochromatic 
single-photon field with an arbitrary spatial shape. If the 
photon is made to pass through a double slit [as shown in 
Fig. 1(a)], its quantum state right after it will be given by 
a 2 x 2 density matrix, p, which can be pure or mixed de- 
pending on the coherence properties of the source and the 
relevant dimensions of the double slit [4, 5, 8-11]. Typi- 
cally, p is written in the basis {|1}, |2}}, which represents 
the slit where the photon was transmitted, so that 
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where Tr(p) = 1, p = p\ and {^\p\4>) > 0, V |V>- 

Taguchi et al. [11] showed that when the photon prop- 
agates through a lens system and is detected by varying 
the transverse position (x) of a pointlike detector between 
its focal plane (far field) and image plane (near field), one 
can access the whole Bloch sphere surface of the spatial 
qubit [Fig. 1(b)]. Figure 1(a) shows a schematic represen- 
tation of this method. A lens of focal length / is placed 
at a distance 2/ from the double slit. The pointlikc de- 
tector can be displaced in the transverse and longitudinal 
directions, x and z, respectively, where x £ (— oo,+oo), 
in principle, and z G [/, 2/]. The effect of a measurement 
over the state p in Eq. (1) for a given position (x, z) is to 
postselect it in a non-normalized state given by 
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\v{x,z)) =^2<pj(x,z)\j), (2) 

3=1 

which, after normalization, corresponds to a point on the 
Bloch sphere surface characterized by spherical polar co- 
ordinates 9 and <j>, as shown in Fig. 1(b). The wave func- 
tion associated with each slit ( j = 1, 2) is given by 

<p s (x, z) = J^e irexd ^'- 1 )/ a sinc {k [x + d(j - A)??]} , (3) 

where A = (D + l)/2, d is the center-to-center separation 
between the slits, a is the slit half-width and sinc(e) = 
(sine)/e. Further parameters are k = 2?:a/Xr]Z, rj = 
(z-f)/f, and Z = (2/ - z)/ V [11]. 

By introducing the operator S(x, z) = 
\tp(x, z))(ip(x, z)\, the detection probability density 
for a photon in the state (1) at the position (x, z) 
becomes (under the assumption that the detector has 
unit quantum efficiency) "P S p a (x, z) = Tt[pS(x,z)]. This 
quantity does not give, directly, information about 
observable quantities for the spatial qubits. In order 
to measure observables here, one must determine the 



FIG. 1: (Color online) (a) Schematic representation of the 
method to measure spatial qubits reported in [11]. (b) Bloch 
sphere: The great circle (in red) connecting |1), |2), and \c±) 
represents the only measurements which do not require com- 
pensation (see text for details). 



longitudinal plane z and two transverse positions xg 
(£ = 1, 2) such that the vectors in Eq. (2) are orthogonal, 
and are thereby the eigenstates of the observable. Then, 
the normalized count rates for these two positions give 
the probabilities for the two possible outcomes labeled 
by I = 1,2 [8-11]. Mathematically, this probability is 
written as 
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PgL = E PWufit' ( 4 ) 

i,3=l 



where <p jt = ipj(xi,z)/JYf i=1 \v>i(x£,z)\ 2 and P s ( pl + 

( 2) 

-Pspa = 1- As an example, at the image plane [z = 2/), 
the Pauli operator a z is measured by detecting in the 
positions xg = (—l) e d/2, corresponding to |^)(^|. 

As mentioned before, this method can be difficult and 
time consuming since it requires the displacement and 
accurate positioning of the detection system for each 
measurement. In addition, when a pair of basis vectors 
is measured in asymmetric positions around the optical 
axis, the detection efficiency between them is different 
due to the diffraction envelope. One must then compen- 
sate for this by using the theoretically expected ratios 
between the corresponding values of the diffraction en- 
velope at those positions. The measurements of observ- 
ables whose eigenvectors lie in the great circle connecting 
|1), |2) and \c±) [red (vertical) circle in Fig. 1(b)] are the 
only ones which do not require compensation. 



III. ALTERNATIVE MEASUREMENT 
STRATEGY 

Here we present an alternative strategy for measuring 
observables for spatial qubits, which is divided into two 
steps. 

A. The POVM and its physical implementation 

The first step of our strategy is the implementation of 
a two-outcome POVM whose goal is to transform, prob- 
abilistically, the initial state of the spatial qubit into a 
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convenient form which will be discussed in the second 
step. The physical implementation of a POVM requires 
the extension of the Hilbert space of the system to be 
measured. This can be provided by an ancillary quan- 
tum system, or ancilla [14]. Here, the ancilla for the 
spatial qubit will be the photon polarization. Therefore, 
the Hilbert space of this larger system is % S pa <£> Upoi, 
where 'H spa (^poi) is the Hilbert space of the spatial (po- 
larization) qubit. Assuming that the systems are initially 
independent, we start by preparing the ancilla in a known 
state, say with horizontal polarization \H). Then, we let 
them evolve subjected to a controlled interaction given 
by a unitary operator U acting on H spa (£> "Hpoi- This in- 
teraction entangles the spatial and polarization degrees 
of freedom and, after it, we measure the ancilla in the ba- 
sis {\H), \ V)} (where \V) denotes vertical polarization). 
The POVM will then emerge as the residual effect on 
the spatial qubit due to the entanglement created by U. 
Finally, we evaluate the postmeasurement state of the 
spatial qubit and the probability for each of the two pos- 
sible outcomes. 

Let lip be the elements of the POVM just described, 
where the subscript p = H,V labels the two possible 
outcomes. Mathematically, they can be written as 



Hp — ApAp 

= (H\U^\p)(p\U\H). 



(5) 



One can easily check that all properties of a POVM are 
satisfied: hermiticity, IP* = tl p ; positivity, (^>|IIj,|^) > 

0, V and completeness, J2 P n P = ^ spa , where I spa 
represents the identity on % S pa- For our purposes the 
unitary operator will be defined as 



U = J2e i <t>i\j}(j\®R(9 j ), 
i=i 



where 



cos(20j) -sin(20j-) 
sin(2^) cos(20j) 



(6) 



(7) 



The operator R(0j), acting on H po i, is written in the 
{\H),\V)} basis. It is easy to see that UU^ = U^U = I, 
where I is the identity on "H spa 8> "H po i • This unitary op- 
eration rotates the photon polarization by 29 j and adds a 
phase shift (f>j, conditional to the passage of the photon 
through the slit j. Therefore, after U, spatial and po- 
larization degrees of freedom become entangled and the 
projective measurement \p)(p\ on the ancilla polarization 
accomplishes the POVM on the spatial qubit. 

If the state of the spatial qubit prior to the mea- 
surement is p and the outcome of the measurement 
is p, then the postmeasurement state of the system is 
p p = A p pA p /Tr(Ilpp). The denominator gives the prob- 
ability {Pp) for the outcome p, or equivalcntly, the prob- 
ability that the state p p has been prepared. Explicitly, 



single 
photon 



1 


1HWPI9,) Pol 


0") 




\ C 




HWP(9 2 ) [j 




U r 

Lens sllt 
@x= 



Det 



FIG. 2: (Color online) Schematic representation of the 
method to measure spatial qubits proposed in this work. 
HWP: half-wave plate; PS: phase shifter; Pol: polarizer that 
transmits ff-polarized light. 



if p = H, the postmeasurement state using Eqs. (1) and 
(5)-(7), will be given by (omitting the polarization) 



Ph = 



1 2 

— J2 Pii****-** cos(2^) cos(2fl i )|i)01, (8) 



and it is prepared with probability 
P H = Tr(n ff/ 5). 



(9) 



An experimentally feasible way of realizing this POVM 
is sketched in the box shown in Fig. 2. First, the photon 
incident on the double slit is prepared with polarization 
\H), so that after transmission its quantum state will be 
p (g> \H){H\ with p given by Eq. (1). Second, a rotating 
half-wave plate (HWP) behind each slit of the double 
slit, followed by a phase shifter (PS) in one of the slits, 
implements the conditional unitary operation given by 
Eqs. (6) and (7), which couples spatial and polarization 
degrees of freedom. Finally, the polarization is measured 
in the basis {\H), \V)} with a polarizer [15]. Clearly, to 
build such a device in a homemade fashion is difficult in 
practice given the very small dimensions (a and d) of the 
apertures (usually, tens of micrometers) . However it can 
be easily implemented with the help of programmable 
liquid-crystal displays (LCDs). LCDs are multipixel op- 
tical devices which can, in a controlled way, rotate the 
light polarization according to the pixel's variable phase 
rctardance. Recently, they have been used to manipulate 
the amplitude [16] and phase [12, 13] of spatial qudits 
with high accuracy, which ensures the feasibility of the 
POVM described here, even for high-dimensional qudits 
as we discuss later. 



B. Photon detection in a fixed transverse position 

The second step of our strategy comprises the photon 
detection in the position (0, z) [17] if the outcome of the 
previous measurement wasp = H (see Fig. 2). Therefore, 
it will be realized with probability Pu given by Eq. (9), 
which makes the whole process probabilistic. In this case, 
the spatial qubit state, which has been transformed into 
Ph [Eq. (8)], is postselected in the non- normalized state 
1) + |2), according to Eqs. (2) and (3). The total prob- 



4 



ability density for detecting the photon will be 
PtotM = P H V spa (0,z) 

2 

= |^(0, z)\ 2 Pije 4 *' coB(20i) 

x e-**' cos(2^), (10) 

where y>i(0, z) = ^(O, z) and is given by Eq. (3). 

By looking at Eqs. (4) and (10) it is easy to check that 
if we set the HWP angles and the phase shift obeying the 
conditions 

fa - fa = fa = arg ( — ] , (11) 

the sum in those equations will be identical. Here I labels 
a possible set of angles (0g,fa). Hence, by choosing two 
sets (£ = 1,2) corresponding to two orthogonal spatial 
projectors, we will be able to reproduce the statistics of 
the observable associated with these basis vectors. De- 
noting the probability for each outcome as P tot , we get 

Aot = Papa) 

such that + P^l = 1. Experimentally, this will be 
achieved by recording the count rate for each set (0£, <pg) 
and then normalizing the obtained data. 

The role of the POVM {n p } (5) in the context of mea- 
surement of observables is now clear. It transforms the 
initial state (1) into state (8) with the angles given by 
Eq. (11), such that its new matrix elements are (up to 
normalization Pg) identical to those elements in the sum 
of the probability (4) for the corresponding spatial post- 
selection. This is done with probability Pg. After that, 
the detection in the point (0, z) enables the sum in the de- 
tection probability density to be independent of <fj(x, z), 
as can be seen in Eq. (10). Therefore, after normalizing 
the count rates we obtain the relation (12). 

Although probabilistic, our method requires neither 
the displacement of the detection system nor compensa- 
tion in the measured count rates, since the detector will 
always be kept fixed at the maximum of the diffraction 
envelope in the point (0,z) [17]. In addition, the possi- 
bility of using programmable LCDs can result in a con- 
siderable reduction in experiment time, since instead of 



positioning the detector for each measurement, we have 
only to update the predefined configuration (fit, fa) in the 
LCD and record the counts for each measurement. 



C. Extension to high-dimensional qudits 

The extension of this method to spatial qudits is 
straightforward. First, assume that in all the previous 
equations where a sum appears, it runs from 1 to D. 
Then, in the unitary transformation (6) we set the polar- 
ization rotations and phase shifters behind each slit such 
that 

e^cos(2^) = ^, Vj,£=l,...,D (13) 

where j labels the slit and I the measurement. After mea- 
suring, in the same way as before, D sets (0ji,fai) cor- 
responding to D orthogonal spatial projectors and nor- 
malizing the data we achieve the relation (12). 

IV. CONCLUSION 

In summary, we have presented an alternative mea- 
surement strategy for spatial qubits which can, proba- 
bilistically, reproduce the statistics of any observable for 
this system. Unlike spatial postselection [11], our method 
requires neither moving the detection system nor com- 
pensation in the count rates. In addition, its extension 
for measuring high-dimensional qudits is straightforward. 
The more challenging step of the strategy would comprise 
the realization of the POVM, but, as suggested here, this 
could be achieved with the help of programmable mul- 
tipixcl liquid-crystal displays. We anticipate that this 
strategy can be a valuable tool in a variety of optical 
implementations of quantum information protocols with 
spatial qubits and qudits such as, for instance, remote 
state preparation, quantum state discrimination, and en- 
tanglement concentration. This will be investigated in 
future works. 
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